(N 









> 



S 



Anisotropic Bianchi-I Cosmological Model in String Cosmology with Variable 

Deceleration Parameter 



Chanchal Chawla^*, R. K. Mishra^ ''' , Anirudh Pradhan^ 

-'^'^Department of Mathematics, Sant Longowal Institute of Engineering and Technology, Longowal— 148106, 

Punjab, India 



^Department of Mathematics, Hindu Post-graduate College, Zamania-232 331, Ghazipur, India 
E-mail: pradhan@iucaa.ernet.in; pradhan.anirudh@gmail.com 



^D ■ *E-mail: c.chawlal37@gmail.com 

(N ■ , . . . 

1,^ ' 'E-mail: ravkmishra@yahoo.co.in 

o 

S: 

Ph. Abstract 

S ■ The present study deals with spatially homogeneous and anisotropic Bianchi type I cosmological model 

representing massive strings. The energy-momentum tensor, as formulated by Letelier (Phys. Rev. D 28: 
2414, 1983) has been used to construct massive string cosmological model for which we assume that the 

^ • expansion scalar in the model is proportional to one of the components of shear tensor. The Einstein's 

field equations have been solved by considering time dependent deceleration parameter which renders the 

^ . scale factor a = (sinh(at))n, where a and n are constants. It has been detected that, for n > 1, the 

presented model universe exhibits phase transition from early decelerated phase to accelerating phase at 
present epoch while for < n < 1, this describes purely accelerating universe which is consistent with 
recent astronomical and astrophysical observations. Moreover, some physical and geometric properties of 
the model along with physical acceptability of the solutions have been also discussed in detail. 
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,3 ! 1 Introduction 

<N . . . . 

The origin of galaxies and formation of other large-scale structures of the Universe are based upon the 
symmetry breaking phase transitions after the big bang explosion, when the temperature falls below some 
critical temperature as predicted by Grand Unified Theories (GUT) (Zel'dovich et al. 1975; Kibble 1976, 
^% . 1980; Everett 1981; Vilenkin 1981, 1985). Such type of transitions lead to the formation of topologically 

stable defects such as vacuum domain walls, cosmic strings, monopoles, textures and other 'hybrid' crea- 
tures. To study the effects of cosmic strings is one of the most outstanding problems in cosmology as the 
vacuum strings (Zel'dovich 1980) generate density perturbations which are strong enough for the formation 
of galaxies. Moreover, the large scale network of strings during the early universe does not contradict the 
present-day observations (Kibble, 1976). These strings also produce gravitational effects that arise from the 
coupling between their stress-energy and gravitational field. The general relativistic treatment of strings was 
pioneered by (Letetier 1979, 1983) and (Stachel 1980). (Letelier 1979) formulated the energy- momentum 
tensor for classical massive strings while (Letelier 1983) represented some cosmological solutions of massive 
strings in Bianchi type I and Kantowski-Sachs space-times. The cloud of strings are termed as massive strings 
and each massive string is formed by a geometrical string with particles attached along its extension. There- 
fore, the cloud forming strings are the generalization of Takabayasi's relativistic model of strings (known as 
p-strings) where particles and strings are together. Since strings are not observable today, so we can say 
that the strings shrink and disappear with the evolution of universe and these cosmic strings end up with a 
cloud of particles. Recently, (Pradhan et al. 2007) and (Yadav et al. 2009) prevailed inhomogeneous string 
cosmological models formed by geometrical strings and used these models as a source of gravitational fields. 
In recent years, cosmic strings have been studied by several authors (Banerjee et al. 1990; Krori et al. 1990; 



Roy and Banerjee 1995; Wang 2005, 2006; Bali and Anjali 2006; Bali et al. 2007; Bali and Pradhan 2007; 
Yadav et al. 2007a, 2007b; Reddy et al. 2007; Pradhan et al. 2008, 2009, 2010, 2011b; Pradhan and Mathur 
2008; Rao et al. 2009, 2011; Pradhan 2007, 2009; Behnchon 2009a, 2009b; Rao and Vinutha 2010; Pradhan 
and Chouhan 2011; Yadav 2012) in Bianchi type space-times in different physical context. 

The simplest of anisotropic models, which, nevertheless, rather completely describe the anisotropic effects, 
are Bianchi type I (BI) spatially homogeneous models whose spatial sections are flat but the expansion or 
contraction rate is directional dependent. The advantages of these anisotropic models are that they have a 
significant role in the description of evolution of early phase of the Universe and they help in finding more 
general cosmological models than the isotropic FRW models. Observations by the Differential Microwave 
Radiometers on NASA's Cosmic Background Explorer recorded anisotropy in various angle scales. It is 
speculated that these anisotropies hide in their hearts the entire history of the cosmic evolution down to re- 
combination, and they are conceived to be indicative of the Universe geometry and the matter composing the 
Universe. The theoretical argument (Misner 1968) and the modern experimental data support the existence 
of an anisotropic phase, which turns into an isotropic one. The isotropy of the present-day universe makes 
the BI model a prime candidate for studying the possible effects of an anisotropy in the early universe on 
modern-day data observations. 

In literature, it is common to use a constant deceleration parameter (Akarsu and Kilinc 2010a, 2010b; 
Amirhashchi et al. 2011a; Pradhan et al. 2011a) as it suggests a power law for metric function or corre- 
sponding quantity. But recent observations from Type la Supernovae (SNe la) (Riess et al. 1998, 2004; 
Perlmutter et al. 1999; Tonry et al. 2003; Clocchiatti et al. 2006) and CMB anisotropies (de Bernardis et 
al. 2000; Hanany et al. 2000; Bennett et al. 2003) revealed that the current universe is not only expanding 
but also accelerating. Several researchers (Pradhan et al. 2006; Pradhan and Otarod 2006; Pradhan et al. 
2012, Yadav 2012) have discussed evolution of the Universe with variable mean deceleration parameter (DP). 
Moreover, (Pradhan et al. 2012) proposed the law of variation of scale factor as increasing function of time in 
Bianchi type VIq space-time, which generates a time dependent DP. This law provides explicit form of scale 
factors governing the Bianchi type VIq universe and facilitates to describe accelerating phase of the Universe. 

Motivated from the studies outlined above, in this paper, the Einstein's field equations have been solved 
for massive strings with time dependent DP in Bianchi type I space-time. This scenario facilitates to describe 
the transition of universe from early decelerated phase to the present accelerating phase. The paper has the 
following structure. The metric and field equations are presented in Section 2. In Section 3, exact solutions to 
the field equations with cloud of strings have been presented with various physical and geometric properties 
of the model. Section 4 proves physical acceptability of the derived solutions and final conclusions are 
summarized in the last Section 5. 

2 The Metric and Field Equations 

We consider totally anisotropic Bianchi type I line element, given by 

ds^ ^-df + A'^dx^ +B^dy^ + C^dz^, (1) 

where the metric potentials A, B and C arc functions of t only. This ensures that the model is spatially 
homogeneous. 

The Einstein's field equations (in gravitational units 8ttG = c = 1) read as 

in-lRgi-T^, (2) 

where T/ is the energy-momentum tensor for a cloud of massive strings and perfect fluid distribution, given 

by 

Ti = {p + p)uiU^ + pgf - XxiX^ . (3) 

Here p and p are, respectively, the energy density and the isotropic pressure for a cloud of strings with parti- 
cles attached to them; A is the string tension density; u' — (0, 0, 0, 1) is the four velocity of the particles and 
X* is a unit space-like vector representing the direction of string. The vectors m* and x* satisfy the conditions 



UiU^ = —XiX'^ = —\,u^Xi = 0. Choosing x* parallel to d/dx, we have a;* — {A ^,0,0,0). 

If particle density of the configuration is denoted by pp, then 

p = Pp + A. (4) 

The Einstein's field equations ([2]) for the line element ^ and energy distribution ([3]), lead to the following 
set of independent differential equations: 

B C BC 

B + c''BC=-P + '^ (5) 

C A cA 

C + A + CA = -P^ (^) 

A B Ab 

A + B + AB=-P' (^) 

Ab bc cA 

AB + BC + CA=P- (^) 

The energy conservation equation T.^ = 0, leads to the following expression: 

p + 3ip + p)H-Xj=0, (9) 

with H being the mean Hubble parameter, which for Bianchi type I space-time can be defined as 

where Hi ~ ^, H2 = ■§ and H^ ~ ^ are directional Hubble factors in the directions of x-, y- and z-axes 
respectively. Here an over dot denotes derivative with respect to the cosmic time t. Also, a is the average 
scale factor of Bianchi type I model written as 

a={ABC)i. (11) 

The usual definitions of the dynamical scalars such as the expansion scalar (6) and the shear scalar (ct) are 
considered to be 

= < = - (12) 



and 



where 
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\fi\ (13) 



The anisotropy parameter (A,„) is defined as 

•^...4i:(^)^ <") 

i— 1 ^ ^ 

The field equations ([5]) — ([8]) involve six unknown variables viz. A, B, C, p, p and A. In order to solve field 
equations explicitly, we require two additional constraints relating these variables, which we shall consider in 
the following section. 



3 Solutions of the Field Equations 

Firstly, we define the mean deceleration parameter (q) as 

ad (h + H^\ ,,, 

^=--^ = -y-^^j=m say. (15) 

The motivation to choose such time dependent DP is behind the fact that the Universe exhibits phase transi- 
tion from the past decelerating expansion to the recent accelerating one, as revealed by the recent observations 
of SNe la (Riess et al. 1998, 2004; Perlmutter et al. 1999; Tonry et al. 2003; Clocchiatti et al. 2006) and 
CMB anisotropies (de Bernardis et al. 2000; Hanany et al. 2000; Bennett et al. 2003). In their preliminary 
analysis, it was found that the SNe data favor recent acceleration (z < 0.5) and past deceleration {z > 0.5). 
From the cosmological observations in the literature (Cunha and Lima 2008; Cunha 2009; Li et al. 2011; 
Frieman et al. 2008; Melchiorri et al. 2007; Ishida et al. 2008; Pandolfi 2009; Lima et al. 2010), the transition 
redshift (zj) from decelerating to accelerating expansion is given by 0.3 < zt < 0.8. More recently, the High-Z 
Supernova Search (HZSNS) team have obtained zt = 0.46 ± 0.13 at (1 a) c.l. (Riess et al. 2004) which has 
been further improved to zt = 0.43 ± 0.07 at (1 a) c.l. (Riess at al. 2007). The Supernova Legacy Survey 
(SNLS) (Astier et al. 2006), as well as the one recently compiled by (Davis et al. 2007), yield zt ~ 0.6 (1 <t) 
in better agreement with the flat ACDM model {zt — (2ilA/^m)^ — 1 ^ 0.66). Now for a universe which was 
decelerating in past and accelerating at the present time, the DP must show signature flipping (Riess et al. 
2001; Padmanabhan and Roychowdhury 2003; Amendola 2003). At present, one may examine the variation 
of cosmic acceleration, i.e., deceleration parameter (DP) with time instead of knowing only the beginning of 
cosmic acceleration and the present value of q. Thus, our choice of variable q is physically acceptable. 

Equation (|15p may be rewritten as 



-+6^ = 0. (16) 

a a^ 

In order to solve the Eq. (J16p . we assume b = 6(a). It is important to note here that one can assume 
b = b{t) = b{a(t)), as a is also a time dependent function. 

The general solution of Eq. p^ with the assumption b = 6(a), is obtained as 

e^^'^'^da^t + k, (17) 

where k is an integrating constant. One cannot solve ()17|) in general as 6 is variable. So, in order to solve 
the problem completely, we have to choose J -da in such a manner that (|17|) be integrable without any loss 
of generality. Hence we consider 

-da = In /(a), (18) 

a 

which does not affect the nature of generality of solution. Hence from ^17} and (fT5|) . we obtain 

f{a)da^t + k. (19) 



Of course the choice of /(a) in (|19p is quite arbitrary, but for the sake of physically viable models of the 
Universe consistent with observations, we here consider 

r7a"~i 

/(«) - ?° .„ - (20) 

ay 1 + a^" 

where a is an arbitrary constant and n is a positive constant. In this case, on integrating Eq. (|19p and 
neglecting the integration constant k, we obtain the exact solution as 

a(i) = (sinh(at))". (21) 

This relation 1^1^ generalizes the value of scale factor obtained by (Pradhan et al. 2012) in connection with 
the study of dark energy models in Bianchi type VIq space-time and the one suggested by (Amirhashchi et 



al. 2011b) in order to study the evolution of dark energy models in FRW universe filled with a mixture of 
barotropic fluid and dark energy by considering time dependent q. 

Secondly, we assume that the component cr^i of the shear tensor {cr-'^) is proportional to the expansion 
scalar (6*), i.e., (t\ oc 9. This condition leads to 



1 /3A B C 



P 



ABC 

I^B^C 



(22) 



which yields 



A B C 



(23) 



_ l+3ff 



where m = ^(i-n\ with /3 being the constant of proportionality. Above equation, after integration, reduces 
to 



A = ki{BCY 



(24) 



where fci is an integrating constant. Here, for simplicity and without any loss of generality, we assume fci = 1. 
Hence we have 

yl=(BC)™. (25) 

The motivation behind assuming this condition is explained with reference to (Thorne 1967), the observations 
of the velocity- red- shift relation for extragalactic sources suggest that the Hubble expansion of the Universe 
is isotropic today within « 30 percent (Kantowski and Sachs 1966; Kristian and Sachs 1966). To put more 
precisely, redshift studies place the limit -jj < 0.3 in the neighborhood of our galaxy today. (Collins et al. 
1980) have also pointed out that for spatially homogeneous metric, the normal congruence to the homoge- 
neous expansion satisfies the condition that ^ is constant. 

Subtracting ([5]) from ([7]), and taking integral of the resulting equation two times, we get 

B 



C 



= ci exp 



C2 [{ABCy^dt 



where ci and C2 are constants of integration. 

Solving ([Tl]) . ([25|) and (|26|). we obtain the metric functions as 

A{t) ==a^, 

B{t) = v^a2('"+i) exp 



'-^iXdt 



C{t) = ^— a2("+i) exp 
^/c^ 



^^ / \dt 



Using the relation ([?!]) in Eqs. (P7|) — ((2^ . we obtain following expressions for the scale factors: 

A{t) = (sinh(Q!i)) "("" 1) , 

C2(-l) 



B{t) = ^/cr(sinh(ai)) 2"(™+i) exp 



1 3 

C{t) = ^=(sinh(ai))2"("+i) exp 



2a 



c2{-iy 



■ cosh(at)F(i) 



■ cosh(ai)i^(t) 



where 



2a 



F{t) = 1 + 1 M + - j cosh^{at) + A M + - j (l + - ) cosh^(ai) + 0[cosh(at)]^ 



(26) 



(27) 
(28) 

(29) 

(30) 
(31) 

(32) 
(33) 



From Eqs. ([50]) — ([5^ . it is clearly observable that the spatial scale factors are zero at the initial epoch t = 0, 
indicating that the model has a point type singularity (MacCallum 1971) with the following form: 



ds' 



'dt^ + (sinh(ai)) "("+!) dx^ + ci(sinh(ai))" 



-u^;-¥T) 



exp 



C2(-l) 2>^ 



cosh(at)F(i) 



dy^ 



1 

Cl 



(sinh(Q!t)) "("+!) exp 



C2(-l) 



3(>» + l) 



■ cosh(at)F(i) 



dz' 



(34) 



For the derived model, the expressions for observational physical quantities such as spatial volume scale 
factor iV), directional Hubble parameters {Hi), Hubble parameter [H), expansion scalar [9), shear scalar 
(u) and mean anisotropy parameter {Am) are given by 



1/ = ABC = (sinh(ai))' 



H.= 



H. 



3a 



n{m + 1) 



coth(at), 



2n{m + 1) 
3q! 



coth(ai) 



C2 



2(sinh(at))* 

C2 



, coth(at) 

2n(m + l) 2(sinh(at)) 

3a 
9 = 3H = — coth(at), 



A^ 



a{2m- 1) 
2n{rn + 1) 



2m- 1 



3th(at) j 



C2 



2(sinh(ai))i 
nc2 



(35) 



(36) 

(37) 
(38) 

(39) 



771+ 1 J 6 Vacoth(ai)(sinh(at))t 

From the above Eqs. (|55 |) — ([55)1 . it is observed that at f = 0, the spatial volume scale factor {V) vanishes 
while the parameters such as scalar of expansion (0), mean Hubble parameter {H) and shear scalar {a) are 
infinite, which is big bang scenario. As t — > oo, V diverges to oo whereas 9, H and a approach to zero. 
Moreover, from Eqs. (1571) and (155)) we obtain 



27n- 1 

2(777+ 1) 



77C2 tanh(at) 
6a(sinh(ai)) n 



(40) 



Clearly 



-> 



i ( 2(m+i) ) as i — ?► OO indicating that ^ -.^ as i — > oo as long as ?77 7^ i. Therefore, the 



dynamics of mean anisotropy parameter {Am) depends upon the value of 777 and for ni 
i — > 00 i.e. the model approaches to isotropy at late times. 



i A 

2 : ^^rJ 



as 



The mean deceleration parameter {q) for the model is found to be 

q = n {1 — tanh {at)) — 1. 



(41) 



For the present age of the Universe i.e. at t — to, q ~ qo, we get the following relationship between the 
constants 77 and a: 



a — — tanh 
to 



1- 



Qo 



1 



(42) 



The sign of q characterizes inflation of the Universe. A positive sign of q i.e. q > 0, correspond to decelerating 
model whereas negative sign of g (particularly — 1 < 9 < 0) indicates accelerating phase or inflationary model. 
From Eq. (|4T|), we observe that q > ior t < itanh"^(l - i)5 and g < for i > itanh^^(l - ^)i . Fig. 
1 depicts the behavior of q with time (i) for different value of 77 and a in accordance with the relation (H^ 
for to = 13.8 GYr and go = —0.73 (Cunha and Lima 2008). It is clearly observable from the figure that for 
< 77 < 1, the presented model universe is in accelerating phase whereas for 77 > 1, the model is evolving 
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Figure 1: The plot of deceleration parameter (q) vs time (i). 

from decelerating phase to accelerating phase. Also, recent observations of SNc la (Riess et al. 1998, 2004; 
Perlmutter et al. 1999; Tonry et al. 2003; Clocchiatti et al. 2006) and CMB anisotropics (de Bernardis et 
al. 2000; Hanany et al. 2000; Bennett et al. 2003) expose that the present universe is accelerating and the 
value of DP lies to some place in the range —l<q<0. 
Using the values of scale factors from (|30 l) ~ ([32|) into the field equations dSJ — ([8]), we obtain the expressions 




Figure 2: The plot of isotropic pressure {p) and energy density (p) vs time {t). Here a = 0.1325, n = 2.75, 
m = 0.6, C2 = 0.00001. 

for string tension density (A), energy density (p) and pressure (p) for the derived model (j34]) as 

3a'^{l-2m) 



^ = V ^n [(3 - ") coth2(ai) + n] , 
2n^(m + 1) 



4to + 1 /3a , , , , 

— coth(ai) I — 



C2 



4n2(m + 1)2 



{m + lf\2n ' 7 V2(sinh(at))^ 

coth^(Q!i) + — -) — ^ coth^(at) + 1 + 



2n{rn +1) V '^ 





(44) 


« yi 


(45) 


2(sinh(at))iy' 



From Eq. (J44|) . we observe that energy density (p) of the fluid is a decreasing function of time and p > 
when 

T^yTT^>^^^i^, (46) 

Say'Am + 1 




Figure 3: The plot of string tension density (|A|) vs time (t). Here a = 0.1325, n — 2.75, m = 0.6. 




Figure 4: The plot of particle energy density (pp) vs time (i). Here a — 0.1325, n = 2.75, m = 0.6, 
C2 = 0.00001. 



where T = sinh(ai). 

Fig. 2 depicts the decreasing behavior of isotropic pressure (p) and energy density (p) versus time (t). Here 
we observe that p is a positive decreasing function of time and it approaches to zero as i — > oo. Moreover, 
pressure (p) varies from early positive values to negative values at late times indicating that the positive val- 
ues of p resembles for decelerating phase of the Universe while the negative pressure indicates some repulsive 
force, which may be responsible for the present accelerating universe which is corroborated by the results 
from recent SNe la observations. 

With the help of expressions (|43l) and (|44|) for the string tension density (A) and the energy density (p), Eq. 
Q yields the value of particle density (pp) as 



Pp 



3a2(8TO2 + 16TO-l) ,n, , 3a2(2TO-l) 
^ ■^■coth^(at) ' ^ ' 



1-71 



coth^(at) + 1 



C2 



2(sinh(ai))i 



(47) 



4n2(TO + l)2 2n{m + 1) 

From the above expression, we notice that pp is a decreasing function of time and it is always positive when 

C2n{m + 1) 



T" 



[3(4m2 + 6m - 1)(1 + T^) - 2n{m + 1)(2to - 1)] ''■ 



> 



\/3a 



(48) 



Fig. 3 and Fig. 4 depicts the decreasing behavior of density parameters (|A| and pp) versus time (t). It is 
clearly observable from the figures that the particle energy density (pp) is positive throughout the evolution 




Figure 5: The plot of f^") vs time (t). Here a = 0.1325, n = 2.75, m = 0.6, C2 = 0.00001. 



m 



lAi 



of the Universe. Moreover, Fig. 5 gives the plot of — versus t. It is interesting to notice that ■'— i- < 1 i.e. 

Pp Pp 

Pp > \X\ during the cosmic expansion (Kibble 1976; Krori et al. 1990), especially during the early universe 
indicating that the early universe is dominated by massive strings (Letelier 1983) evolving with deceleration 
but at later times it got disappear. Since there is no direct evidence of strings at present epoch, therefore, we 
are interested in the construction of only that models of the Universe which evolves from the era dominated 
by geometric strings or massive strings, followed by particle dominated era with or without the remnants 
of strings. Therefore, the derived model describes evolution of the Universe consistent with the present day 
observations. 



It is also observed that the above set of solutions satisfy the energy conservation equation © identically. 
Therefore, the above solutions are exact solutions to the Einstein's field equations ([5]) — ([8]). 

4 Physical Validation of the Model 
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Figure 6: The plot of deceleration parameter (q) vs redshift (z). Here n = 2.75 and qo ~ —0.73. 
Here, we firstly define the relationship between average scale factor (a) and redshift parameter (z) as 



-1 + ^, 



(49) 



where oq is the present value of the scale factor i.e. at z = 0. 



Using the value of scale factor from Eq. (|2ip into the above relationship, one can easily obtain the 
expressions for cosmic time (t), mean Hubble parameter (H) and mean deceleration parameter (q) in terms 
of redshift parameter (z) as follows: 



t(z) = — sinh 
a 



j n-l-qo 
((Zo + l)(z + l)2"' 



(50) 



H(z) = - coth 
n 



sinh 



1 



qo 



{qo + l){z + l)^ 



q{z) = n — 1 — n 



-I 2 



tanh sinh 



n-1 



90 



(go + l)(^ + l) 



2n 



(51) 



(52) 



Such type of relation ([5 
expansion for q[z) = go 



provides a two-parameter (n,qo) parametrization just like a linear two-parameter 
■ qiz (Riess et al. 2004), where go is the present value of DP and qi is the deviation 
in the redshift evaluated at z = 0. If we set go = —0.73 (Cunha and Lima 2008), then a positive Zt may be 
obtained only for the positive values of gi since go is negative and the dynamic transition (from deceleration 
to acceleration) occurs at q{zt) = 0, or equivalently, Zt = — — . Another parametrization of considerable 
interest is q{z) = go + gi2;(l + z)^"^ (Xu and Liu 2008) where the parameter gi describes the total correction 
in the distant past {z ^ 0, q{z) — go + 9i) and positive zt may be obtained for the positive values of gi > |go|. 
Also, the dynamic transition occurs at Zt = — „'^°.„ ■ Moreover, (Xu et al. 2007) have also proposed one 
such parametrization g(z) = 1/2 — a/(l + z)^ where a and h are constants which can be determined from 
the current observational constraints. As z — >■ oo, g — > 1/2 which resembles the value of DP at dark matter 
dominated epoch. In this case, the present value of DP is 1/2 — a and zt = (2a)^/^ — 1. 

For the derived model, the expression for transition redshift {zt) is found to be 

1 
n - go - 1 



Zt 



-1 



,(n-l)(go + l). 

Fig. 6 shows the variation of g with redshift {z) for go — —0.73 and n 
Zt = 0.351817 which is in fine tuning with (Xu et al. 2007). 



(53) 



2.75. It clearly elaborates that 



4.1 Physical Acceptability of the Solutions 




18 20 



Figure 7: The plot of energy conditions vs time t. Here a = 0.1325, n ~ 2.75, m = 0.6, C2 = 0.00001. 

Firstly, we check for the fulfillment of energy conditions for the derived model. The weak energy condi- 
tions (WEC) and dominant energy conditions (DEC) are given by 



10 



{i) p>0, 

{ii) p + p> 0, 

{Hi) p — p> 0. 

The strong energy condition (SEC) is given hy p + 3p> 0. 

The left hand side of energy conditions have been depicted in Fig. 2 and Fig. 7. From both the figures, 
we observe that 

• The WEC and DEC are satisfied in our model throughout the evolution of universe. 

• The SEC is satisfied during the early stages of the evolution of universe whereas it violates at present 
epoch. The violation of SEC gives anti-gravitational effect due to which the Universe gets jerk and the 
transition from the earlier decelerated phase to the present accelerating phase takes place (Caldwell et 
al. 2006). Hence the presented model clearly describes the late time acceleration of the Universe. 



p ' 0,8275 




Figure 8: The Plot of sound speed (vg) vs. time t. Here a = 0.1325, n — 2.75, ni — 0.6, C2 = 0.00001. 

Secondly, we check for the physical acceptability: velocity of sound (vg) should be less than the velocity of 
light (c). As we are working in the gravitational units with unit speed of light, i.e., c = 1, therefore the 
velocity of sound {vg) must satisfy the condition < Vs{— ■^) < 1. 

From Eqs. (H^ and P5|) . we obtained the sound speed as 

dp o? [2n(TO + l)(2m + 1) - 3(4771^ + 2to + 1)] (sinh(ai)) ""^ ^ „c|(m + 1)^ 



dp 



3a2(4m + l)(sinh(ai))^"2 _ nc?^{m + 1)2 



(54) 



Fig. 8 depicts the plot of sound speed (w^) with time (i). It is clear form the figure that the speed of sound 
is almost constant throughout the evolution of universe and Vg <\. 

On the basis of above discussions and analysis, we conclude that the derived solutions are physically 
acceptable. We also emphasize here that the behavior of physical parameters depicted in the figures are due 
to a selected choice of model parameters while it may change for some different sets of parameters. 

5 Concluding Remarks 

In this paper, exact solutions to Einstein's field equations for spatially homogeneous and anisotropic Bianchi 
type I model with string fiuid as a source of matter and variable DP have been established. The choice 
of variable DP yields average scale factor a(t) — (sinh(Q!t))n which represents two types of models: the 



11 



accelerating models (0 < n < 1) and the models with transition (n > 1) from the early decelerated phase to 
present accelerating phase (see, Fig. 1) which is in good agreement with the recent supernovae observations. 
For the presented model (for n = 2.75), the Universe undergoes a phase transition from positive pressure 
(deceleration) to negative pressure (acceleration). We conclude that this negative pressure plays the role of 
repulsive force which may be responsible for the current accelerating phase of the Universe. It is also observed 
that the early universe is dominated by massive strings evolving with deceleration but eventually disappear 
from the Universe at late time i.e. at present epoch. Moreover, the derived solutions are physically acceptable 
in concordance with the fulfilment of WEC and DEC and Ug < 1 while the violation of SEC is reproducible 
with current astrophysical observations. Finally, we conclude that the model represents shearing, non-rotating 
and expanding universe, which starts with a big bang and exhibits point type singularity at i = 0. Thus, 
the solutions demonstrated in this paper may be useful for better understanding of the characteristic of 
anisotropic Bianchi type I model in the evolution of the universe. 
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